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Abstract. We validate a method of computing spin Chern numbers that was developed 
with Hastings and only conjectured to be correct. Specifically, the Pfaffian Bott index can 
be computed by the "log method" for commutator norms up to a specific constant. 

We establish lower bounds on the distance almost commuting unitaries must be moved 
to change an index. This is done for the Bott index, the winding number index for general 
unitary matrices, and for the Pfaffian Bott index for self-dual matrices. 
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1. Introduction 



<^ ' We begin with a survey, and some improved theorems, in the subject of almost commuting 

unitary matrices and the winding number index discovered with Exel [6] . Our eventual goal is 
to improve known mathematical results regarding invariants of almost commuting matrices 
by making them quantitative, so that they be more useful in applications, especially in 
relation to topological insulators. In the present paper we stick with stability results, with 
the expectation that in numerical studies of index the averaging over disorder [151 II2J E] 
can be done more intelligently. However, we anticipate applications in the classification of 
\ C*-algebras. 

The seasoned veterans from the world of C*-algebras are warned that the Bott index 
discussed here is the one developed with Hastings to be computable numerically, and it 
involves different functions in its definition from what is usual [61 Definition 1.2]. That such 
invariants, contingencies related to stable relations, are not canonically defined is a known 
difficulty, discussed in [3]. 

We follow mathematical conventions, so U* refers to the conjugate-transpose when U is a 
matrix. 

Our main result, from the standpoint of phyics, is that the spin Chern number of a fi- 
nite system in 2D can be computed as follows. As discussed in length in [12] . the position 
operators, assume the geometry of the model is the torus, are two commuting unitary ma- 
trices, but the band-compressed position operators are almost commuting matrices acting 
on low-energy space that are almost unitary. If the diagonalization of the Hamimtonian is 
done preserving symmetry, if we start in class All we end up with matrices A and B so that 
A" = A and B$ = B (within machine precision) and only have approximate relation 

AB « BA 

A* A ps / 

B*B ps / 

If we compute any approximate self-dual logarithm iK of B, meaning 

— 7T < K < 7T 
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and K} = K and e lK B, then the spin Chern number of the system will be even or odd 
depending on the sign of 

Where {•, •} is the anti-commutator and Q is a specific matrix discussed below that creates 
anti-symmetry in the formula so that that Pfaffian makes sense. 

Of course the logarithm, as a matrix function, will be discontinuous as we cannot assume 
anything about the spectum of A and B beyond Kramers doubling. Put another way, the 
various self-dual approximate logarithms of B are far apart, so it is not even clear that the 
sign of the Pfaffian is well defined. This issue we adrress in the following sections. 

There is a separate issue of how to compute approximate logarithms of almost unitary 
matrices, and how to be sure to get a self-dual output given a self-dual input. That is 
discussed in a separate paper [T3] . 

In practice, it is feasable to simply compute the unitary parts of A and B before computing 
further, so 



U = A(A*A)~* 

which is unitar and self-dual, as is 

V = B(B*B)~*. 

The comparison of||[{7, V]|| with \\[A, B]\\ is not hard. For the rest of the paper we discuss 
self-dual matrices. We discuss the Pfaffian-Bott index, which depends only on U and V (or 
A and B) but was shown [121 Lemma 5.8] to equal the spin Chern index, at least for large 
systems and other mild assumptions. 

Our results on the Pfaffian-Bott index are quantitative, but the correspondance with spin 
Chern numbers is not quantitative. This will be possible once certaining new quantita- 
tive results concerning localization of Wannier functions or concerning approximation by 
commuting matrices are developed. 

Our results are principally stated in terms of unitary matrices. However, the study of 
almost commuting unitary element of C*-algebras is not that different. We know this because 
we know that the soft-torus is RFD [2]. This needs to be checked for the real soft torus, 
whose definition is easy enough, now that we have a theory of generator and relations for 
real C*-algebras [20] . 

A unitary U in a C*-algebra A might have a logarithm. Given K* = K in A so that 
U = e tK and ||[?7,V]|| < | for a second unitary element V in A, then find that the Bott 
index in K (A) will be determined by 



That is, Theorem 16.21 holds true for all unital C*-algebras. Indeed, the proof is the same as 
in the matrix case. 
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2. Phantom and ordinary almost commuting matrices 
Given two unitary matrices U and V with 5 = \\ [U, V] ||, we find 

\\VUV*U* - I\\ = \\[U,V]\\ 

and so by the spectral theorem 

a(VUV*U*) C {z e T\\z- 1| < 5} . 

Thus when 5 < 2 we can define {VUV*U*) t for t between and 1, using a branch of x l with 
discontinuity on the negative a>axis. This is a continuous path of unitary matrices from I 
to VUV*U* and 

det(VUV*U*) = det(J) = 1 

so 

t ^ det ((UVITV*)*) 

is a loop on the unit circle. We define w{U, V) to be the winding number of this path. 

This winding number invariant is very computable. There are a few alternate formulas, 
including 



(2.1) 



to (U, V) = Tr ( — log (VUV*U*) 



which we easily prove: since in some basis VUV*U* is diagonal and unitary, 

vuv*u* 



for some — tt < 6j < n and 



Tr ( — log (VUV*U*] 



i9 1 



2ni 



-Tr 



and 



Jt6i 



det {{VUV*U*f) = det 
and this path has winding number 



o it0„ 



Lemma 2.1. W^/ien [/ anc? V are commuting unitary matrices, ou(U, V) = 0. 
Proof. In this case the path of determinants is the constant path. 

Lemma 2.2. For 

( 1 \ / e ^/« 



1 
1 



■. 

1 0/ 



V 



-in/n 



□ 
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we have u(U, V) = — 1. 
Proof. We findVW*l7* = e~^r I and so 

1 1 / 27T? 

TWlog(lW*[/*)) = Tr / ) = -1. 



□ 



It is easy to modify the example in Lemma 12.21 to get a pair of unitary matrices with 
|| [U, V] || =5 and u(U, V) = n for any < 5 < 2 and any n. 

Theorem 2.3. Pairs of unitary matrices with \\[U,V}\\ = 5 and lu(U,V) ^ exist for all 
< 5 < 2. For swc/i a pa«r, the distance to a commuting pair of unitary matrices exceeds 1, 
meaning 

\\u-u4 + ||v- Ki|| > i 

whenever U\ and V\ are unitary matrices with U\Vi = V\U\. Indeed, 



l^-t/JI + llV-Vill > \ l-\5 2 



The proof of this will be broken into lemmas and propositions. Theorem 12.31 is a slight 
improvement on the results obtained with Exel [6]. 

Lemma 2.4. Suppose \\[U 0l V Q ]\\ < 2. If \\[Ui, Vi\\\ < 2 and 

u(U Q ,Vo) ^ uiUuVx) 
then for any path U S ,V S of unitary matrices, there must be at least one so so that 

\\[U S0 ,V S0 }\\=2. 

Proof. We will use a homotopy argument. If no such So existed then 

(s,t)^det((u s v s u:v;f) 

is a homotopy between the path that determines co(U ,Vo) and the path that determines 
u(Ui, V\). Therefore the winding numbers of these paths are equal. □ 

Proposition 2.5. Suppose \\{U,V]\\ < 2. If ||[tfi, Vi]|| = 2 or if \\ [U u V x ] \\ < 2 and 

then 

\\u-UtW + ||v-Vi|| > yfi-\\\[u,v]\\ 2 + yfi-iWiWf. 

Proof. The general case follows from the case ||[E/i,Vi]|| = 2 by Lemma 12.41 Since this 
invariant in zero for commuting pairs, this proposition implies Theorem 12.31 

Since —1 is in the spectrum of V\U{V*Ul then, by a spectral variation result pQ, similar 
to Weyl's estimate on Hermitian operators, 

WVxUiVfUl - VUV*U*\\ > dist (-l,a(VUV*U*)) . 

We also have 

WViUxVfU* - VUV*U*\\ < 2 \\U - U x \\ + 2 \\V - V x \\ 

and 

5 = dist(l,a(VUV*U*)) . 
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Figure 2.1. The solid curve shows the minimum radius in the gap at —1 in 
with spectrum of VUV*U* for unitaries U and V, as a function of 5 = \\[U,V}\\. 
The lower curve show the minimum distance one must go to find a pair with 
winding number index either undefined or different. 

A little basic geometry shows 

(dist (-1, a {VUV*U*))f + (dist (-1, a (VUV*U*))) 2 = 4 

and so 

\\U-Ux\\ + \\V-Vx\\ > ^l-\5 2 . 

□ 

We now get to a difficult question. Is the winding number invariant the only obstruction 
to closely approximating U and V by commuting unitary matrices? It is important here that 
we stick with the operator norm in defining "close approximation" as the answers to these 
sort of questions can change dramatically if considering the Frobenius norm [8j [9j [HI EH] • 
(In particular, see the discussion in section III in [JS]-) Results such as this also change 
dramatically when the matrices come from different symmetry classes, as seen in [T6] . 

There is an answer, but it is only a non-quantitative, nonconstructive result for small 6. 
This is Theorem 6.15 in [3] and we repeat it here. Also it matters that we are only interested 
in results for unitaries in M rf (C) that are independent of d [TOT [13]. 

Theorem 2.6. For any e > 0, there is a 5 in (0, 2) so that, whenever U and V are unitary 
matrices in M^(C) with \\[U,V]\\ < 5 and u(U,V) = 0, there exist unitary matrices U\ and 
Vi in M d (C) so that 

\\U-UxW + \\V -V x \\ < e 

and {U U V X \ = 0. 

A serious limitation of the invariant u(U, V) is that the only known proofs of Theorem 12.61 
work with another invariant k(U, V) that is only equal to u(U, V) for small 5. This invariant 
has become known as the Bott index. 
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Another limitation of u(U,V) is that we don't know how to modify it to work in other 
symmetry classes. For example if we have self-dual unitary matrices, so Z7" = U and = V, 
where (j is a specific generalized involution detailed below, we find 

(VUV*U*f = u*v*uv 

and so generally VUV*U* is not self-dual. 

3. A DIRECT fT-THEORY INVARIANT — THE BOTT INDEX 

We need functional calculus of unitary matrices, also called matrix functions in applied 
mathematics. An example is above where we applied the logarithm to a unitary matrix. 
Generally speaking, for the functional calculus /(V) to be define for a unitary matrix we 
need / defined on the circle. One diagonalizes V via another unitary Q and applies / on the 
diagonal, so 

/ (e iei ) 

V = Q I ••. I Q* =► f(V) = Q I •-. I Q* 



e i8d 



However, most of our calculations will involve Fourier series, and traditionally those are 
defined in terms of scalar functions that are periodic. 

Definition 3.1. Assume then that / is periodic of period 2ir we define /[V] as /(V) where 
/ = / o log. In other words, 

/ \ /fie,) 

V = Q\ ••. ]Q* f[V\=Q\ •-. \Q*. 

\ e"'/ V f(0i. 

When / has uniformly convergent Fourier series, this is easier: 

oo oo 

(3.i) f{x) = a ^ nx =► f[v] = a -y n - 



Definition 3.2. Define 



and g and h: 



and 



f(x) = — — (150 sin(x) + 25 sin(3x) + 3 sin(5x)) 
128 



'o ^e[-|,f] 



h 



V^P [-f,|] 

which are shown in Figure 13.11 For any unitaries set 



muv) -f f\V\ g[V] + l{h[V],U}\ 

1 U ' V} ~ \g\y\ + l{h\y],u*} -f[v] ) 



QUANTITATIVE ^-THEORY AND SPIN CHERN NUMBERS 



7 




-3-2-10123 




-3-2-10123 -3-2-10123 



FIGURE 3.1. Functions for the standard Bott index (trig method). 

and if \\[U, V}\\ < 0.206007 define k(U, V) as the integer 

K(U,V) = ^Sig(B(U,V)). 

For an invertible, hermitian matrix A we define its signature Sig(A) as the number (with 
multiplicity) of positive eigenvalues minus the number of negative eigenvalues. We need to 
know that \\[U, V]\\ < 0.206007 forces B(U, V) to be invertible. 

It should be noted that 

7i = f(0 2 ) 

l2 = g{9 2 ) + h{9 2 )cos{e x ) 
73 = /i(0 2 )sin(0i) 

defines the coordinates of a map from T 2 — > S 2 C M 3 that has mapping degree one. Also 
notice gh = and f 2 + g 2 + h 2 = 1. 

Theorem 3.3. Suppose U and V are unitaries and 

S = \\[U,V]\\ < 0.206007. 

(1) The hermitian matrix B(U,V) has a spectral gap at of radius at least ^|a/1 — 55. 
Indeed, the gap is at least as large as the function of 5 plotted as a solid curve in 
Figure HO . 

(2) The distance \\U — Ui\\ + || V — Vi|| needed so that B(U\, V\) has in its spectrum is 
at least |yl — 55. Indeed, this distance is at least as large as the function of 5 plotted 
as a dashed curve in Figure HO . 
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Figure 3.2. In Both: Solid curve on top is the guaranteed spectral gap in 
B(U, V) for a given 5. The dashed curve is the distance one can move where it 
is proven the gap will not close. The dotted curve in the left plot is ^jVl — 55. 
The dotted curve on the right is |a/1 — 55. 

This will be a lot of work to prove, and the gap is much smaller than we saw for VUV*U*. 
Why do we bother? The point is symmetry. 

Suppose U* = U and V" = V for unitaries in M2d(C) and with fl the generalized involution 
discussed in the next section, that physicists call the dual. Then B(U,V) is in M 2 d(C) £g> 
M 2 (C) which has on it the generalized involution r = In terms of real C*-algebras, this 
is a copy of M M+2 (C) with the transpose operation. In physics language, we are tensoring 
two half-odd-integer spin systems to get a system with integer spin, in a non-standard basis. 
We find 

B{U,V) T = -B(U,V) 

and so B(U, V) defines a class in 

K 2 (R) = tf_ 2 (H) = Z/2 

that is computed directly in terms of the a Pfaffian, hence the Pfaffian-Bott index studied 
in [15]. In return for a small gap, indeed no guaranteed gap if || [U, V] || is too large, we get a 
construction that is amenable to symmetries. 

An easy upper bound on the gap radius can be found, using the example in Lemma I2.2[ 
shown in Figure [331 This shows we cannot get a big a gap using B(U, V) as was possible 
with VUV*U*, but that the situation is likely not as bad as Figure [3721 indicates. 

Fortunately B(U,V) is readily computable since /, g and h we chosen to have rather 
fast decay in their Fourier coefficients. Thus we were able to replace (13.11) by the simpler 
evaluation of order-5 trig polynomials. For applications to index studies, the following is the 
most useful. We will later have a version of this for the Pfaffian-Bott index. 

Proposition 3.4. Suppose \\[U,V}\\ < 0.206007. // ||[L7i, 14)11 < 0.206007 and 

k{U u Vr) ^k{U,V) 
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Figure 3.3. Upper bound on the gap radius, as ||[V, U]\\ varies from to 
about 0.21. 

then 

Wu-UtW + wv- v.w > ^i-5\\[u,v]f + ^i-sii^^ii 2 . 

Another consequence of Theorem 13.31 is the following. 
Theorem 3.5. Suppose U and V are unitary matrices. If || [U, V}\\ < 0.206007 then 

k{U,V) = u(U,V). 

Proof. Exel [5] showed that for ||[C/,V]|| < 2 the winding invariant equals an abstract K- 
theory invariant. In the notation of [5J, this is defined in terms of b e in K (A e ). It is easy to 
check that as long B(u, v) has a spectral gap, where u and v are the generators of the soft 
torus A e , the i^-theory class of B(u, v) is b e . □ 



4. Pfaffians 



The Pfaffian of skew-symmetric matrices is not the most familiar object, and it it not 
clear at the outset how it applies to a problem involving self-dual matrices. Let us start by 
recalling the dual operation. 

We fix 

/ 



-I 



in M 2 at(C) and this specifices the dual operation 



ZX T Z 



as 



" A 


B ' 


8 




-B T 


C 


D 






A T 
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When we discuss M 2 (M 2A r(C)) = M 2A r(C) <g> M 2 (C) = we require the unitary 



Q 



which has the convienient property (T2J Lemma 1.3] 









/ 








-il 


1 


/ -iZ ' 


1 





I 


il 





71 


iZ I 


"71 





il 


I 











-il 








/ 



Q*X m Q = (Q*XQ) 



Here 



or 



" A 


B ' 


8<X>B 






C 


D 




-0 





A n 


A 12 


A 13 


A u 


m 


A T 


— ^34 


-A T 

^24 


A T 


A 2 i 


A 22 


A 23 


A 24 




-A T 

^43 


A T 

^33 


^23 


-A T 


Aai 


A 32 


A33 


A 34 




-A T 

^42 


^32 


^22 


-A T 


A 41 


A 42 


^43 


A M 






^31 


-A T 


A T 



Recall the Pfaffian is defined for all skew-symmetric, complex 2n-by-2n matrices by 



/ 



Pf 



O 



\ 





— CLi 



a 1 


-02 



«2 



-03 



03 





\ 



J 



det(0)aia 3 • ■ ■ a 2n -i- 



/ 



for O real orthogonal. (All skew-symmetric matrices have such a factorization, a modified 
Hessenberg decomposition.) The essential properites are that 

Pf(YXY T ) = det(Y)Pf(X) 

for arbitrary Y, that the Pfaffian varies continuously, and 

(Pf (X)) 2 = det(X) 

so is zero exactly on the set of skew-symmetric, singular matrices. 

For matrices with the symmetry X"®" = —X we can define a modified Pfaffian 

Pf(X) = Pf (Q*XQ). 

We still have 

(pf(X)) 2 = det(X) 

and that this varies continuously. The sign of the Pfaffian can be used to prove a homotopy 
result, in the same way we use the determinant to detect that the real orthogonal matrices 
fall into two connected parts. 

Proposition 4.1. Suppose B is in M 4 tv(C) and B* = B and B T = — B and B is invertible. 
Then 

Pf(B) Gl\{0}. 

If B\ and B 2 are elements of 

n= {B e M m (C)\ B* = B = -B T is invertible} , 
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then they can be connected by a path in % if and only if with Pf(Si) and Pf(_B 2 ) have the 
same sign. 



Proof. We can apply Theorem 8.7 in [12J to iB and learn that there is a real orthogonal 
matrix O so that det(O) = 1 and 



B = 



zAi 
-i'Ai 



i\ 2 
-i\ 2 



6 ] 



and all the real numbers Aj are positive except Ai which as the same sign as Pf(5). It is 
clear from this form that two such matrices with Pfafhan of the same sign will be connected. 
The spectrum of B is {±Ai, . . . , ±A2tv}. Thus there will be an even number of negative 
eigenvalues, so det(-B) > 0. Since the square of the Pfaffain is the determinant, we find 
Pf(.B) is real, and as B is invertible, the Pfaffian cannot be zero. Since the Pfaffian varies 
continuously, it is not possible to connect two matrices in H that have Pfaffians of opposite 
signs. □ 

Now we explain the Pfaffian-Bott index. 

Definition 4.2. Let /, g, h and B(U,V)be as in g3j If \\[U,V}\\ < 0.206007 define k(U, V) 
as the value in {±1} given by 



k 2 (U,V) = Sign (Pf(B(U,V))^ 
Lemma 4.3. When U and V are commuting unitary matrices, w(U, V) = \. 



Proof. One proof is to use Lemma 6.4 of jTT]. Here is another. 

One easily checks that n 2 remains constant along a path so long as \\[U t , V t ]\\ < 0.206007. 
One can use functional calculus, and so keep the self-dual condition, to deform a commuting 
pair Uq and Vq that is self-dual over to XJ\ = I and V\ = I. One can chen compute that 



and 



Pf 



/ 
-/ 



Pf 



B(I,I) 



iZ 
-iZ 



/ 
/ 



Pf (iZ) Pf (-iZ) 



□ 



Proposition 4.4. Suppose \\[U, V]\\ < 0.206007 and that U ,V , U\, V\ are self-dual unitary 
matrices. If \\[U u Vi]\\ < 0.206007 and 



K 2 (U 1 ,V 1 ) ^ K 2 (U,V) 



then 



u - Ux\\ + \\v - k|| > o v 1 - 8 v \\\ + o v 1 - 8 R> Vi 
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Proposition 4.5. Suppose \\[U,V}\\ < 0.206007 and that U,V, XJ\, V\ art self-dual unitary 
matrices. If Ui commutes with V\ and k 2 (U,V) = — 1 then 

\\U-U 1 \\ + \\V-V 1 \\>± + yi-8\\[U,V]\\ 2 . 

We get larger, if more complicated, estimates below. 

5. Proofs 

We are able to estimate the radius of the spectal gap in B(U, V) by esitmating the norm 
of B(U,V) 2 -I. 

Lemma 5.1. Suppose f, g and h are continuous, real-valued functions that are 2n -periodic, 
and with 

f 2 + g 2 + h 2 = l 

and 

gh = 0. 

Suppose U and V are unitary matrices and define 

f[V] g[V] + \{h[V],U} 

g[V] + ±{h[V],U*} -f[V] 



Then S* = S and 



\S 2 -l\\<2\\[h[V},U}\\ + \\[f[V],U}\\ 



Proof. Since /, g and h are real- valued, the matrices f[V], g[V] and h[V] are hermitian. Let 
us write / for f[V], etc. We see easily S* — S and 



/ 

g+HKU*} 



9+\{h,U} 
-f 



1 2 


" I 


" 




' A 


B ' 







I 




B* 


A* 



where 



and 



We have 



A = f 2 + g 2 -I + \{h,U}{h,U*} + lg{h,U*} + l{h,U}g 
= -h 2 + \ {h, U} {h, U*} + \g {h, + \ {h, U} g 

B = fg+\f{h,U}-gf-\{h,U}f 
= \f{KU}-\{h,U}f. 



S 2 -I < 



A 

A* 

+ \\B\ 



+ 



B 

B* 



Notice f 2 + g 2 + h 2 = l forces these functions to take value in [—1, 1] so ||/[V] || < 1, etc. 
Therefore 

< \ \\hUhU* + Uh 2 U* + UhU*h - 3h 2 \\ + \ \gU*h - ghU* + hUg - Uhg\\ 

< \ \\h\\ \\Uh - hU\\ + \ \Uh 2 - h 2 U\ + \\g\\ \\XJh - hU\ 

< 2 \\Uh - hU\\ 
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and 

\\B\\ = ±\\hfU + fUh-hUf-Ufh\\ 
<l\\h[f,U]\\+l\\[LU}h\\ 
<\\[fM\ 

so 

\\S 2 - I\\ < 2 || [h,U] || + || [f,U] || . 

□ 

Definition 5.2. Suppose / is continuous and 27r-periodic. Following [17] we define rjf : 
[0,oo)^[0,oo)by 

Vf (S) = sup {|| [f[V], A] || | V unitary, ||A|| < 1, \\[V,A]\\ < 5} 

where the supremum is taken over all V and A in every unital C*-algebra. 

Once we have a bound on rjf we can use it to bound more than just commutators. Indeed, 
by [T7] Lemma 1.2], for any two unitaries V^and V2 we have 

\\f\y]-m]\\<vt (\\v -v.w). 

We need a special case of a Lemma in |17j . 

Lemma 5.3. Suppose f, g andh are continuous and periodic, and that g is the trigonometric 
polynomial 

n 

9( x ) = Yl bkeikx - 

k=—n 

then rj(S) < m5 + b where 

n 

m = \kbk\ 

k=—n 

and 

b = max h(x) — min h(x) 

We need this only in the case of real-valued functions, in which case b is the diameter of 
(the range of) h. 
A computer algebra package is recommended here, as it shows 

407 / 96 9 \ 

fix) 2 H cos 6 (x) 1 H cos (2x) H cos Ux) = 1 

JK J 512 K J \ 407 v ; 407 V ' ) 

which means 

, . /407 „ , / 96 " " 9 ~~ / . \ 

g{x) = V 5T2 C ° S V 407 C ° S + 407 C ° S V ~ ^ ) 

s /407 3 . . r 96 9 ~~ , . 

K x ) = V 512 C ° S 407°° + ^ cos ( 4x )^[-f-f]( x ) 

A handy formula here is 

407 / 96 , n . 9 tA A / 15 ... 9 4 . A 

1 + — — cos (2x) H cos [Ax) = 1 H cos (x) + — cos (x) 

320 V 407 v ; 407 v ' J \ 40 v ' 40 v 7 
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and we get alternate expression for g and h, in particular 



K x ) = cos 3 (a;)^l + ^ cos 2 (x) + cos 4 (x)x[-f, §](>)• 



We new bound the derivative of (7 and /i, computing 

To 



15 



9 



cos (x) \ 1 + — cos 2 (x) + — cos 4 (a; 



40 



40 



p(sin(x)) 
16-\/g(sin(x)) 



where 



and 



p(x) = 30x(x - l)(x + l)(3x 4 - 10x 2 + 15) 



q(x) = 9x 4 - 33x 2 + 64. 
On [—1, 1] the max of p(x) is 150 and the min of q(x) is 64 so we find 

150 



\h'(x)\ < 



128 



and the same for g' . 

We need /iasa cosine or Fourier series so need 



in . . /407 3. , / 96 , n , 9 , . . , 

Cn = ^ / cos(nx) V^ cos (x) V +407 C ° S ^ ' + 407 C ° S ^ 

"2 

We computed these with numerical integration, and without checking error estimates, in 
[15] . We compute these a little more carefully here. We find 



1 n 
2^ 



407 



cost //./• ) \/ — — cos 3 (x) 

fc=0 



1/2 

jfe 



— — cos (2x) H cos (Ax) ) dx 

407 v ' 407 v M 



fc=0 

00 



V— / 407 A/ 2 



2tt V 512 V fe 

5> 



96 



9 



cos(nx) cos (x) I — — cos (2x) + —— cos (4x) dx 



407 



407 



fc=0 



where the I n> k were defined in-line and are easy to compute with a computer algebra package. 
The convergence here is rather rapid, as 



2tt 



1 /407,1/2 



512 V k 



96 



9 



cos(ra) cos (x) I — — cos (2x) H cos (Ax 



407 



407 



dx 



< im(l/2\ (_105\ k 
- 2tt V 512 V jfe y V 407/ 

Letting T# denote the the Taylor polynomial 



T K (x) « Vl + a; 
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15 



n 





1 


2 


3 


4 


5 


a n 





15Ui 

256 





256 





3i 
256 


b n 


0.202047 


-0.179940 


0.125655 


-0.066010 


0.023445 


-0.003886 




0.202047 


0.179940 


0.125655 


0.066010 


0.023445 


0.003886 



Table 1. These are approximations to the first coefficients in the Fourier 
expansions of the /, g and h used to define the Bott index. Extend these to 
negative indices by the rules a_ n = a^" and 6_ n = b n and c n = —c n . 



of degree K expanded at 0, we have 
^ i n,k\ - 27r y 512 



k=K+l 



k=K+l 
oo 



9-ttV K19 I 



2tt V 512 

1 /407 
27V 512 



.fc=0 



1/2 
k 



105 
"407 

105 
407 



K 

E 

fc=0 



1/2 



105 
407 



1 105 T 
1_ 407 " Tk 



105 
407 



This means we need K — 7 to get six digits absolute accuracy: 
Using the values in the table to define 

5 

h 5 (x) = c n e inx 



n=— 5 



or 



hs(x) = Co + 2ci cos(x) + 2c2 cos(2x) + 2c3 cos(3x) + 2c4 cos(4x) + 2cs cos(5x) 
and we find, 



A. 

dx 



(h - h 5 



150 

< + 1.48498 = 2.656855. 

~ 128 



and so we can estimate to six decimal places the maximum of \h — h 5 \ by simply plugging 
in values between — 7r and 7r with an even spacing of a little less than 10 -7 . Keeping track 
of the errors and rounding up, we find 



diamO(x) - h 5 {x)) < 0.004110 



and we note 



5\\F 



E 

n=— 5 



n\c r , 



1.48498. 



The other estimates of this sort, for ho, ■ ■ ■ , h±, are summarized in Table |2j 
We also can use brute force to find 

sup \h(x) - h 5 (x)\ < 0.002338. 

X 

Lemma 5.4. For any unitary matrix V , 

\\h 5 [V] - h[V}\\ < 0.002338 
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n 


m = bound on \\h' n \\ F 


b = bound on diam (h(x) — h n (x)) 








1 


1 


0.359880 


0.732237 


2 


0.862500 


0.350141 


3 


1.258560 


0.106619 


4 


1.446120 


0.017509 


5 


1.48498 


0.004110 


oo 


2.99208 






Table 2. Bounds on rjh as a slope and an offset. 



We get the same error estimate on using only 6_ 5 through 65 when numerically computing 

g[v). 

We also use the folklore estimate 



where 



\[hlV\,U\\\<\\h%UU,V\\\ 

11*11, =E 



\nc r , 



is in terms of Banach algebra norm || • || F on functions whose Fourier series converge absolutely, 
where \\q\\ F is the £ 1 -norm of the Fourier series of q. 

Lemma 5.5. For h as in Definition 13.21 we have 



Proof. We check that 



\h% < 2.99208. 



|== (45 cos 4 (x) + 60 cos 2 (x) + 120) 



Z3256 



h'(x) = 

^1 + ^cos (23;) + ^ cos (Ax) 
and attack this as three factors. It is easy to see 



sm(z)cos (xjX[-f,f]W 



1 



V3256 

and the next factor is not so bad, as we see 



(45 cos 4 (a;) + 60 cos 2 (2) + 120) 



225 
V3256 



1 



1 + H cos (2x) + cos (Ax 



00 

k=0 
00 

=E 

k=0 



-1/2 
k 

-1/2 
k 



96 . , 9 

cos (2a;) H cos (Ax) 

407 v ' 407 v ' 



105 



(_!)« _ 
v ; 1 407 



1 _ 105 
407 



407 
302' 



We estimate 



sm(x)cos 2 (x)x[-i,%](x) 
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n 


m = bound on \\f' n \\ F 


b = bound on diam (f(x) — f n {x)) 








2 


1 


1.171875 


0.21875 


2 


1.7578125 


0.0234375 


oo 


1.875 






Table 3. Bounds on rjf as a slope and an offset. 



as follows. The Fourier series of — isin(x) cos 2 (x)x[-z-,z-] is 

5 w 5 o -1 r ' ' -1 nr ? -1 ^ ? -1 r ' -1 ' ' 1 n ? -1 r ' -1 r> ' -1 nr ? ? n-i r ' ? 



3465vr 315vr 105vr 16 15vr 16 16 15vr 16 105tt 315tt 3465vr 
with terms beyond n = 3 being given by 

/7rn\ / 1 1 

~ C ° S \2~) V(^-l) 3 -4(n-l) + (n + l) 3 - 4(n + 1) 

Therefore 

sin(x) cos 2 (x)x[ 



" 2 ' 2 J 



I 16 18 4 A 1 



1 1U IO <± ^ — V 

"4 15tt 105tt (2n + l) a - 4 (2n + 1) 

1 16 18 4 16 4 f°° 1 

< — I 1 1 1 h - / — 5 = dx 

~ 4 15vr 105vr 315vr 3465vr vr J 4 8x 3 + 12x 2 - 2x - 3 

1 16 18 4 16 1 , /81 

= — I 1 1 1 1 In — 

4 15vr 105tt 315tt 3465tt 4vr V 77 



and so 



^ 225 /407 /1 16 18 4 16 1, /81\\ 

\\h'\L < \ - H 1 1 1 1 In — < 2.992076 

~ V 302 \4 15tt 105tt 315tt 3465vr 4tt \77 ) ) 

□ 

We approximate / the same way, but this is just arithmetic since / is already a trigono- 
metric polynomial. 

Lemma 5.6. Let f and g and h be as in Definition ^. 2 . Then 

T]f(8) < m5 + b 

for each of the values in Table 0, and 

Tj g {5) < m5 + b 

and 

T]h(S) < m8 + b 

for each of the values in Table 
Let 

which is shown in Figure ISTTl 
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0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



Figure 5.1. The function /3(5) that bounds \\B(U, V) 2 - I\\ in terms of 5 = \\[U, V}\\. 
Theorem 5.7. Suppose U and V are unitary matrices. Then 

\\B{u,vf-i\\<nmv]\\) 

d for V]|| < 0.206007 the gap at in the spectrum of B(U, V) has radius at least 



an 



The other key thing we have is how B(U, V) varies. 
Theorem 5.8. If Uj and Vj are unitary matrices then 

\\B {U , V ) - B (U , V )\\ <P{\\V - Vi\\) + \\U - U x \ 

and so 

\\B (U , V ) - B (U , V )\\ < {\\V - T4|| + \\U - U Y \\ 
All our main theorems now follow. 



6. The log method 

An alternate way to compute the Bott index was considered in [6]. One replaces B(U, V) 
with 



Bl(U,V) 




where iK is the logarithm of V, meaning — ir < K < tt and e lK = V. Numerical evidence 
in [15] suggests that, for small commutators, the Pfaffian-Bott index can be computed using 
B L (U, V). We validate this here. 

Since the logarithm is not continuous, numerical errors will mean we might accidentally 
compute the wrong branch of logarithm on V, or indeed any logarithm at all. 
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We note that when q is periodic, 

q(K) = q[V\. 

Lemma 6.1. Suppose f, g and h are real-valued Borel functions on [— satisfying 

f 2 + g 2 + h 2 = l 

and 

gh = 0. 

Let q(x) = f(x)h(x) and assume further that q and h are continuous and 2n -periodic, 
pose U and V are unitary matrices and —iK is a logarithm of V and define 



S = 



f(K) g(K) + ±{h(K),U} 

g(K) + ±{h(K),U*} -f(K) 



Then S* = S and 



P 2 - I\\ < (\\ 9 \\ + 1) || [Mn + \ W\h[VlU]\\ 2 + \ || [h 2 [V],U] || + \\[q[V],U]\\ . 

Proof. We write / for f(K) and compute a bit more carefully than before. We find 

1 

4 



hUhU* + Uh 2 U* + UhU*h - 3h 2 



and 



= - \\hUhU* -h 2 + UhU*h - Uh 2 U* + 2Uh 2 U* - 2h 2 \ 
= J \\-[h, U][h, U}* + 2 (Uh 2 U* - h 2 
<]\\[h,U]\\ 2 + h\[h 2 ,U}\\ 



l -\\ g {h,U*} + {h,U}g\\= l -\\gU*h + hUg\\ 



\g[U\h\ + [h,U\g\ 



and 



<IMI 11^,^11 



\ \\f{K u} - {h, u}f\\ = ± \\fhu - ufh + fuh - huf\\ 

= \ \\2fhU - 2Ufh + fUh - fhU - hUf + Uhf\\ 

< \\fhu - uhf\\ + \ \\fuh - fhu\\ + \ \\huf - uhf\\ 

< \\fhu - uhf\\ + \ \\uh - w\\ + \ \\w - uh\\ 
= ||[A,f/]|| + || [h,u] || 

= II [?, t>] II + P,t>] II- 
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-3-2-1 1 2 3-3-2-1 1 2 3 




-0.2 - 

-0.4 - 
-0.6 - 
-0.8 - 



-3-2-10123 



Figure 6.1. First part of the path, showing f t then h t then g t . 




-3-2-1 1 2 3-3-2-1 1 2 3 



Figure 6.2. Second part of the path, showing only f t then h t . Here gt is zero. 

Theorem 6.2. Suppose U and V are unitary matrices, i/ || [£/, V] || < | then for any choice 
of K with —7i < K < 7i and e lK = V, there is a path B t of invertible self-adjoint matrices 
between B{U, V) and 




and, ifU and V are self-dual, then the path may be chosen with the symmetry B\ = B t . 
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Proof. We can select paths ft, gt and h t from 

(fo,ga,h ) = (f,g,h), 



the standard triple, as in Definition 13. 2\ to (fi,g%, hi) where 



(6.1) fi(x) = -x, g x {x) = 0, /ii(x) = \/l ^x 2 



B t (U,V) 



-x, gi{x) = 0, h-i(x) = \ 1 - — 

7T V TT 

and the conditions f 2 + g 2 + h 2 = l and gh = hold along the way. This gives us paths of 
matrices 

f t (K) g t (K) + l{h t (K),U} 

g t (K) + \{h t (K),U*} -f T (K) 

with the needed symetries. It remains to show these are invertible. One needs to compute 

(6.2) (N| + 1) + \ (vkM) 2 + \v h ^) + vM 

and check that this takes value less than 1 at 5 = ~. This is a bit much to do by hand, so we 
instead show by computer that (16. 21) takes value less than 0.95 at S — g for all t in a mesh 
t\, . . . , t w selected so that 



\\f t . - ft-., II +|k--0t- +1 || +ll^ -^ +1 || < Vl - 0.95 w 0.2236. 
We can keep / t fixed at 

-1 -7T<X<-| 

/ t (x) = { (150 sin(x) + 25 sin(3x) + 3 sin(5x)) -f < x < f 
1 If < a; < 7T 

while altering ^ from the standard g to 0. The more interesting part of the path interpolates 
/ f from the above to while keeping g t = and 

/i t (x) = v 7 ! - /t(ar). 

The graphs of the computed bounds are shown in Figure 16.31 These bounds have been 
rounded up to accomodate the various errors in computing offset terms when applying 
Lemma 15.31 The errors in computing Fourier coefficients lead to sub-optimal results, but do 
not need to be accounted for as it is the computed coefficients that are used when applying 
Lemma 15.31 The analysis of the error bounds is dull and omitted. 

It is apparant that that the limitation on the constant in this result comes from the 
functions used in the log method (16. ip . The computed bounds are shown in Figure [631 D 

The computations of the Pfaffian-Bott index are mostly validated by Theorem 16.21 What 
remains is validaing the algorithms used to compute a self-dual logarithm of a self-dual 
unitary. This is discussed in [Lf] . 
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